Useful multiparticle entanglement and sub shot-noise sensitivity 
in experimental phase estimation 
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We experimentally demonstrate a general criterion to identify entangled states useful for the 
estimation of an unknown phase shift with a sensitivity higher than the shot-noise limit. We show 
how to exploit this entanglement on the examples of a maximum likelihood as well as of a Bayesian 
phase estimation protocol. Using an entangled four-photon state we achieve a phase sensitivity 
clearly beyond the shot-noise limit. Our detailed comparison of methods and quantum states for 
entanglement enhanced metrology reveals the connection between multiparticle entanglement and 
sub shot-noise uncertainty, both in a frequent ist and in a Bayesian phase estimation setting. 



The field of quantum enhanced metrology is attract- 
ing increasing interest [lj and impressive experimen- 
tal progress has been achieved with photons [2] [6], 
cold/thermal atoms [7], ions [8j and Bose-Einstein con- 
densates [9j[T0]. Several experiments have demonstrated 
phase super resolution (3[ [8], which, if observed with a 
high visibility of the interference fringes, allows to utilize 
the state for quantum enhanced metrology [5j[TT]. So far, 
only few experiments have implemented a full phase esti- 
mation protocol beating the shot-noise limit AO = 1/ VN 
with N > 2, where N is the number of particles 0- 
9J. Recently, it has been theoretically shown that sub 
shot-noise (SSN) phase sensitivity requires the presence 
of (multi-) particle entanglement p~2| [T3j. In this letter, 
we experimentally demonstrate this connection. For an 
entangled state and a separable state with N = 4 ad- 
dressable photons, we measure the quantum Fisher in- 
formation (QFI) Fq [14J, which quantifies the amount of 
entanglement of the state useful for SSN interferometry 
|12| . We then show how this entanglement can indeed be 
exploited by implementing a Maximum Likelihood (ML) 
and a Bayesian phase estimation protocol, both clearly 
yielding SSN phase uncertainty. 

The usefulness of an experimental state can be quan- 
tified by the quantum Fisher information (QFI) Fq[p, J] 
[13] . A probe state p of N qubits is entangled and allows 
for SSN phase estimation if the condition 

F Q [p,J]>N (1) 

is fulfilled [12] . Here J = \ J2iLi ls the linear gener- 
ator of the phase shift, and &ty = Hi • a is a Pauli matrix 
rotating the qubit i along the arbitrary direction Hi. The 
maximal Fq further depends on the hierarchical entan- 
glement structure of the probe state and genuine multi- 
particle entanglement is needed to reach the Heisenberg 
limit OUS], the ultimate sensitivity allowed by quan- 
tum mechanics. With N = 4 qubits, 2-particle entangled 



states have Fq < 8, while for 3-particle entangled states 
Fn < 10 [T3] fT6]. The ultimate limit is F Q < N 2 = 16 
which is saturated by the so-called Greenberger-Horne- 
Zeilinger (GHZ) state [H H7J HE] . 

A state fulfilling Eq. ([I]) allows for SSN phase uncer- 
tainty due to the Cramer-Rao theorem, which limits the 
standard deviation AO of unbiased phase estimation as 

A0 > i 1 > / 1 = ■ (2) 

yj mFfi [0 O , p, J] y mF Q [p, J] 

The first inequality defines the Cramer-Rao lower bound 
(CRLB). Here Oo is the true value of the phase shift, m is 
the number of repeated independent measurements, and 

( 3 ) 

The Fisher information Fjx [0q, p, J] depends on the con- 
ditional probabilities P(p\0o) to obtain the result p in a 
measurement when the true phase shift is equal to Oq. It 
is bounded by the QFI [I2j [T3], the equality being sat- 
urated for an optimal measurement /i op t- From Eqs ([T]) 
and ([2| and from the bounds for multi-particle entangle- 
ment, we can infer that, if the experimentally obtained 
Fq of a 7V-qubit state exceeds the value for /c-particle 
entanglement, one can achieve a phase sensitivity better 
than that achievable with any (k — l)-particle entangled 
state of any N qubits [16J. 

For the experimental demonstration, we use the 
symmetric four-photon entangled Dicke state [2TJ [22] 
\Df ] ) = (\HHVV) + \HVHV) + \HVVH) + \VHHV) + 
\VHVH) + \VVHH))/y/6 and the separable state 
l^sep) = | + + + +) observed from multiphoton paramet- 
ric down conversion E3] [Fig. |T]a)]. Here \HHVV) = 
|#}i ® \H} 2 ® \V) 3 ® \V)i, \H)i (\V)i) refer to the hor- 
izontal (vertical) polarization of a photon in the spatial 
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Figure 1: a) Experimental setup. The source uses pulsed 
parametric down conversion with a type II cut /3-Barium- 
Borate crystal (\ pU mp = 390nm) |15) . After passing an inter- 
ference filter (IF), the photons are symmetrically distributed 
into 4 spatial modes by using 3 non-polarizing beam split- 
ters (BS). The Dicke state \D^) is observed if one photon is 
detected in each of the four output arms [23 . The separa- 
ble state iV'sep) is created by inserting a |+) polarizer before 
the first BS. Each polarization qubit is addressed individually 
and rotated by exp[— ia y 0/2] (violet box) by a halfwave-plate 
(HWP). Each polarization analyzer (PA) is composed of a 
quarterwave-plate (QWP), a polarizing beam-splitter (PBS) 
and an avalanche photo-diode (APD). b) Schematic of our 
interferometric setup. 



mode z, and |±) = -^(\H) db \V)). From the measured 

density matrices (p^ p and p^ [15]) we deduce a fidelity 

of 0.8872±0.0055 for \D^) and 0.9859±0.0062 for |^ S e P ) 
(errors deduced with Poissonian count statistics) and also 
the QFI determining the suitability of the experimen- 
tally observed states for phase estimation. For the ideal 

(2) 

Dicke state \D\ ; ), the QFI reaches its maximum value, 
Fq[\D { n /2 \ A = N(N + 2)/2 = 12, when a Hi = & y for 
all z (J = J y ) [24]. In the experiment, this choice leads 
to F Q [p e * p , J y ] = 9.999 ± 0.095, at the maximal value 
achievable with 3-particle entanglement. An optimiza- 
tion over the local directions fii [24], leads to the slightly 
higher value F° pt [p^ p , J°^} = 10.326 db 0.093, detecting 
useful 4-particle entanglement with 3.5 standard devia- 
tions. Sure enough, using a witness operator it is possible 
to prove 4-particle entanglement in a simpler way [T5| f2T] . 
With only a subset of the tomographic data we obtain a 
witness expectation value of —0.2205 db 0.0055, proving 
4-particle entanglement with a significance of 40 stan- 
dard deviations [15J. However, witness operators merely 
recognize entanglement, whereas our criterion directly in- 
dicates the state's applicability for a quantum task. The 
separable state |^sep) ideally allows for sensitivity at the 
shot-noise limit, FQ[\ip sep ), J y ] = N = 4. The experimen- 
tal density matrix leads to F Q [p^, J y ] = 3.894±0.023, a 
value close to the expected separable limit (the optimized 
value being F° pt {p e * p , ^° pt ] = 4 - 014 ± 0.025). 

In order to demonstrate that the precision close to 
the one predicted by Fq can indeed be achieved in 
practise, we experimentally implement a phase estima- 
tion analysis with the input states p e ^ p and Pse P - Our 



interferometric protocol transforms the probe state by 
U(6o) = exp[— i J2t=i using the halfwave-plates 

and phase shifts depicted in Figs [I] a) and b). The un- 
known value of the phase shift 0q is inferred from the 
difference in the number of particles, 2p = Njj — Ny 
(p = —2,-1,0,1,2), in the states V and H. For the 
ideal states and the rotation directions Hi = y, this mea- 
surement is optimal, and hence F^ = Fq. Experimen- 
tally, the optimized direction and measurement can be 
different because of noise and misalignment. However, 
for the observed states the expected improvement would 
be rather small. 

The relation between the phase shift and the possi- 
ble results of a measurement is provided by the con- 
ditional probabilities P(p\0o). These are measured ex- 
perimentally and compared with the theoretical ones for 
both the separable and the entangled state, as shown in 
Fig. [2] a)-k). A fit to the measured conditional prob- 
abilities provides P exp (/i|0), which are used to calculate 
the Fisher information according to Eq. ([3| [see Fig.[2]l)]. 
As expected, our experimental apparatus can surpass the 
shot noise limit for a broad range of phase values (where 
F| xp > 4), and can even exploit useful three particle en- 
tanglement (where i^ xp > 8). 

The phase shift 0q is inferred from the results, 
Pi ,M2, ...,/i TO , of m independent repetitions of the inter- 
ferometric protocol. We will refer to such a collection of 
measurements as a single m-experiment. In the experi- 
ment, we set the phase shift to 9 known values 0$. For 
each #0, 12000 results pi are independently measured and 
grouped into vectors of length m to perform the phase 
estimation for different values of m (= 1, 10, 100). Using 
this data, we implement a ML and a Bayesian phase esti- 
mation protocol. While both have been recently used in 
literature for phase estimation [4j|25], here they are com- 
pared in detail and applied for the first time to demon- 
strate SSN phase uncertainty with more than two parti- 
cles. To display the quantum enhancement and to com- 
pare the methods we use the rescaled uncertainty A res 
defined below. 

In the ML protocol, the estimator est of the unknown 
phase shift is determined as the value maximizing the 
likelihood function C(0) = UT=i P e^i\°) i 19 J- For dif- 
ferent m-experiments it fluctuates with standard devi- 
ation A# est , which has to be calculated by repeating a 
large number of single m-experiments. For large m, the 
distribution of # e st approaches a Gaussian centered on #0 
and of width A# est saturating the CRLB, Eq. Q [19J. 

Fig. [3] shows the distributions of the estimator est for 
the phase shift 0q = 0.2tt and different values of m. As 
expected, with increasing m, the histograms approach a 
Gaussian shape with standard deviation A# est decreasing 
as 1/ ^mFfr. The width of the histograms is smaller for 
the Dicke state (red lines) than for the separable state 
(blue lines). Fig. [4] shows A res = ^/mAQ est as a func- 



3 




Figure 2: Calibration curves and Fisher information. The small panels show the conditional probabilities P e ^ P (p\0o) for the 
state p e £ p [red curve, upper row a)-e)] and for p^ p [blue curve, lower row f)-k)]. Dashed black lines are the ideal probabilities 
P(/jl\0)j dots are experimental results. The red and blue curves are fits obtained by assuming that the main source of errors are 
misalignments in the polarization optics [15] . The measurements are performed for 31 values of Oq by collecting approximately 
7000 events for each phase value. Panel 1) shows the Fisher information [Eq. (|3j], obtained from the fits P e ^p(p\0o). The line 
widths correspond to the error intervals F A [6> ,p exp ] ± AF A with AF^[p e * p ] < 0.08 and AF^pt^} < 0.04 (Tg]. Horizontal lines 
indicate limits for separable states ("sep", equal to ^[l^sep)]), for 2- and 3-particle entangled states, and for the ideal Dicke state 
F/ilJD^)]. Theoretically, Fp, = Fq holds for the ideal input states, phase operations and output measurements. Experimentally, 
we observe F^^o,p exp ] < Fq[p gxp ] due to technical noise. In particular, Fp,[0o, p e ^ p ] is strongly reduced for values of the phase 
shift where some of the ideal conditional probability densities P(p\0o) of panels a)-e) go to (0o = 0, 0.37T, 0.57r). For reduced 
visibilities (when P e ^ P (p\0o) > while ideally P(p\0o) = 0), the contribution to the Fisher information is reduced since in these 
points also the derivatives of P e x P (p\0o) vanish, cf. Eq. pi. 



tion of 0q . For m = 10 the standard deviation is below 
the CRLB (Eq. [2| for several 6$ values. This is possible 
because the estimation is biased, i.e., for b = (0 est ) — 6o 
we have 6^0 and do b ^ [T5l [T§] . The bias can 
be taken into account by replacing the numerator in the 
CRLB Eq. (|| by |1 - d 0Q b\. For even smaller m, only 
few different maxima of the likelihood functions C{6) can 
occur, see Fig. [3] a). Then, est scatters significantly and 
hardly allows for an unbiased phase estimate. When 
m = 100, the bias is strongly reduced and the agree- 
ment of A# est with the unbiased CRLB is improved sig- 
nificantly. While the bias is still large enough to cause 
apparent sensitivities below the shot-noise limit for the 
separable state, for the Dicke state the CRLB is satu- 
rated for a large phase interval. This clearly proves that 
the multiparticle entangled Dicke state created experi- 
mentally indeed achieves the SSN phase uncertainty pre- 
dicted by the CRLB Eq. ^ using the experimentally 
obtained Fisher information from Fig. |2]l). 

A conceptually different phase estimation protocol 
is given by the Bayesian approach assuming that 
the phase shift is a random variable. The proba- 
bility density for the true value of the phase shift 
being equal to 6, conditioned on the measured re- 
sults /ii, /i2, Mm, is provided by Bayes' theorem, 

p(o\{»i}r=i) = Pex P ({»i}r=i\o)p(o)/p({»i}r=i)- to de- 

fine the a priori probability density P{6) we adopt the 
maximum ignorance principle and take P{0) to be con- 
stant in the phase interval considered. The Bayesian 
probability density is then given by i) cx 

ni^i -Pexp(/^|0) = £(0). The phase shift can be es- 



timated as the maximum of the probability density as 
before. However, in contrast to the ML method, the 
Bayesian analysis allows to assign a meaningful uncer- 
tainty to this estimate even for a single m-experiment 
and biased estimators. This can be taken, for instance, 
as a confidence interval AO = C around the estimate, 
where the area oiP{e\{iii}f =1 ) is equal to 68% (see Fig. [3] 
d) and [15]). 

Figs [3] d)-f) illustrate how the Bayesian probability 
density evaluated for a single m-experiment becomes 
Gaussian with a width 1 / y/mi 7 ^, already for small val- 
ues of m. In contrast, the ML histograms [Figs [3] a) -c)] 
approach a Gaussian shape more slowly. We also inves- 
tigated how the Bayesian analysis performs on average 
using the same data as in the ML case. The results are 
shown in Fig. [4] with the rescaled Bayesian uncertainty 
A res = y/mC for various 0q and averaged over several 
m-experiments. For m = 10 the mean values of the con- 
fidences deviate from the CRLB and have a large spread. 
For m = 100, however, the confidences agree well with 
the CRLB for most values of 0q> f° r both states. 

In conclusion, we have investigated experimentally the 
relation between SSN phase estimation and the entangle- 
ment properties of a probe state. We have identified use- 
ful multiparticle entanglement by determining the quan- 
tum Fisher information from the tomographical data of 
a four photon Dicke state. The benefit of such entangle- 
ment has been demonstrated by implementing two dif- 
ferent phase estimation analyses, both of which saturate 
the Cramer Rao bound and clearly surpass the shot noise 
limit. The approach is completely general: it applies for 
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Figure 3: Comparison of the ML method to the Bayesian 
approach for the estimation of a phase shift 0o = 0.2n (ver- 
tical dashed black line). Upper row: histograms (normalized 
to one) of the estimators es t obtained for large number of 
repetitions of ra-experiments: a) m = 1, b) m = 10, and 
c) m = 100. Red (blue) solid lines show the results of the 
state p e ^ p (ptep)- Lower row: exemplary Bayesian proba- 
bility densities P^K/i*}™^) of single m-experiments for d) 
m=l, e) m=10, and f) m=100 for the state p e £ p (solid red 
lines) and ptep (solid blue lines). In panel e) the dashed 
black lines are Gaussians of width 1 / ^/rnF^ plotted to illus- 
trate that the densities rapidly approach a Gaussian shape. 
For m = 1, we plot P(% = 1) for p e * p and P(% = -1) 
for ptep- The shaded region indicates the confidence interval 
[0est — C, Oest + C] around the maximum of the distribution. 



any probe state, is scalable in the number of particles 
and does not require state selection. Our study thus pro- 
vides a guideline for the future technological exploita- 
tion of multiparticle entanglement to outperform current 
metrological limits. 

We thank N. Kiesel, W. Laskowski, and O. Giihne for 
stimulating discussions. R.K., C.S., W.W., and H.W. 
acknowledge the support of the DFG- Cluster of Excel- 
lence MAP, the EU projects QAP and Q-Essence, and 
the DAAD/MNISW exchange program. W.W. and C.S. 
thank QCCC of the Elite Network of Bavaria and P.H. 
thanks the ERC Starting Grant GEDENTQOPT. 



[5] 
[6] 



[9 
[10 

[11 
[12 

[13 

[14] 

[15 
[16 

[17] 
[18 
[19 
[20 



[21 
[22 



[23 

[24] 

[25] 



T. Nagata, eta/., Science 316, 726 (2007). 
Xiang, G. Y., et al, Nature Phot. 5 43 (2010). 
J. Appel eta/., PNAS 106, 10960 (2009); M. H. Schleier- 
Smith eta/., Phys. Rev. Lett. 104 073604 (2010). 
D. Leibfried eta/., Science 304, 1476 (2004). 
C. Gross eta/., Nature 464, 1165 (2010). 
M. F. Riedel eta/., Nature 464, 1170 (2010). 
K.J. Resch eta/., Phys. Rev. Lett. 98, 223601 (2007). 
L. Pezze and A. Smerzi, Phys. Rev. Lett. 102, 100401 
(2009). 
P. 



G. 



Hyllus eta/., |http://arxiv.org/ abs / 1006 .4366 
Toth, |http://arxiv.org/abs/1 006.4368 

C. W. Helstrom, Phys. Lett. 25 A, 101 (1967); S.L. Braun- 
stein and CM. Caves, Phys. Rev. Lett. 72 3439 (1994). 
See supplementary material for additional information. 
A pure state is called /c-particle entangled if \ip k ~ ent ) = 
®z^i wnere \^i) is a non-factorizable state of Ni < k 
particles, and Ni = k for at least one / [15] . 

D. M. Greenberger, M.A. Home and A. Zeilinger, Going 
beyond Bell's theorem (Kluwer Academics, 1989). 

V. Giovannetti, S. Lloyd and L. Maccone, Phys. Rev. 
Lett. 96, 010401 (2006). 

H. Cramer, Mathematical Methods of Statistics (Prince- 
ton Univ. Press, 1946). 

Note that it is possible to further improve the sensitivity 
by applying the phase shift several times to the probe 
system [T8] . 

N. Kiesel eta/., Phys. Rev. Lett. 98, 063604 (2007). 
For indistinguishable qubits the symmetric Dicke state 
reduces to a Twin-Fock state, see M.J. Holland and K. 
Burnett, Phys. Rev. Lett. 71 1355 (1993). 
W. Wieczorek eta/., Phys. Rev. Lett. 103, 020504 (2009). 
P. Hyllus, O. Giihne and A. Smerzi, Phys. Rev. A 82, 
012337 (2010). 

L. Pezze eta/., Phys. Rev. Lett. 99, 223602 (2007); Z. 
Hradil eta/., Phys. Rev. Lett. 76, 4295 (1996). 




SNL 



0.1 0.2 0.3 0.4 0.5 



0.1 0.2 0.3 0.4 0.5 

QJn 



* Present address: Vienna Center for Quantum Science 
and Technology, Faculty of Physics, University of Vienna, 
Boltzmanngasse 5, A- 1090 Vienna, Austria 
^ Present address: Dep. of Theor. Phys., The Univ. of the 
Basque Country, P.O. Box 644, E-48080 Bilbao, Spain 

[1] For a recent review see V. Giovannetti, S. Lloyd, and L. 
Maccone, Nat. Phot. 5, 222 (2011). 

[2] J.G. Rarity eta/., Phys. Rev. Lett. 65, 1348 (1990). 

[3] M.W. Mitchell, J.S. Lundeen and A.M. Steinberg, Nature 
429, 161 (2004); P. Walther eta/., ibid. 429, 158 (2004) 

[4] M. Kacprowicz eta/., Nature Phot. 4 357 (2010). 



Figure 4: Rescaled phase uncertainties A res obtained with 
the probe state p e £ p (red) and ptep (blue), with a) m = 10 
and b) m — 100. The solid red (blue) line is the expected 
uncertainty given by the CRLB [Eq. Q] using the exper- 
imental Ffx [see Fig. ||] 1)] for the state p e * p (ptep)- Stars 
are the results of the ML analysis (standard deviation of 
#est) with A res = ^/m/\0 est - Circles are the results of the 
Bayesian analysis with A res = ^/^(C), and the error bars dis- 
play the scatter of C. Horizontal lines are the shot- noise limit 
(SNL) (solid line), the limit for an ideal Dicke state (D ld ) 

A res = 1/ \J Fq[\D^}] (dashed line) and the Heisenberg limit 
(HL) Ares = 1/4 (dot-dashed line). 



